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Abstract 

The next generation of submillimetre/millimetre in- 
struments will provide us with a deeper insight into the 
mechanisms that rule galaxy formation. As the brightest 
starbursts are thought to be heavily obscured at optical 
wavelengths, the opening of this new window will com- 
plement the present observations, and enable a detailed 
investigation of the hierarchical merging of galaxies at re- 
mote epochs. In this context, we will present preliminary 
results of simulations currently developed for deep galaxy 
surveys with FIRST. 

A'^-body simulations of galaxy mergers are being devel- 
oped to produce realistic morphologies and star formation 
histories. A Schmidt law is used to account for the evolu- 
tion of the stellar and gas content of each particle. These 
simulations will be gathered in a library, including the 
temporal information with a timestep of 10 Myr. We then 
build consistent spectra accounting for this star forma- 
tion history at each resolution element of the simulations. 
These morphologies will be included in the framework of 
the GalICS hybrid model of hierarchical galaxy forma- 
tion, which reproduces the main observational constraints. 
We intend to present synthetic maps with the characteris- 
tics of the FIRST PACS/SPIRE instruments and discuss 
the optimal strategy for deep surveys complemented with 
the ground-based ALMA project. 
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1. Introduction 

The understanding of galaxy formation and evolution is a 
rapidly evolving field. The infrared/submillimetre sources 
detected with SCUBA/JCMT are compatible with the 
IR Background detected with COBE. The dusty galaxies, 
whose identification (redshift, morphology) is so difficult 
(e.g. Downes et al. 1999), probably escape current optical 
surveys. The new generation of instruments, which will be- 
come available in the coming years, will provide us with 
an unprecedented multi-wavelength view of these galaxies. 
From this perspective, we have undertaken a comprehen- 
sive simulation work that will be used to reproduce current 
observations and to predict forthcoming ones. 



In this paper, we present the main framework of this 
study as well as preliminary applications to FIRST. In 
Section we present a simple analytic model of starburst 
spectra and apply it to known starburst galaxies. In Sec- 
tion ^, the simulations of galaxy mergers being developed 
are described together with preliminary applications. In 
Section ^, we conclude with a short overview of the gen- 
eral framework of the GalICS hybrid model of hierarchi- 
cal galaxy formation, within which the previous work will 
be incorporated. 
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Figure 1. Best model for the spectra of the ERO 
JI645O2+46 26.4 (HRIO). T he measurements have been 



published by Dey et al. ( 1999, 



2. A SIMPLE MODEL FOR GALAXY SPECTRA 

We have built an analytic tool to produce spectra of star- 
bursts, that will be used to produce images of galaxy 
mergers based on N-body/SPH simulations. Each of these 
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spectra will be associated with each resolution element, 
and it is essential to have a restricted number of parame- 
ters. This minimal treatment is intended to model high-z 
galaxies and to get a deeper insight into the process of 
star formation. 
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Log A (mm) sensitivity for 



5a, 111 exp. time 



Figure 2. A typical strong starburst galaxy (HRIO) spec- 
tra is extrapolated for different redshifts. The thresholds of 
PACS, SPIRE as well as NGST and ALMA are presented. 
Such galaxies could potentially he observed up to a redshift 
of 5. 



2.1. Assumptions 

First of all, we gathered a library of Simple Stellar Pop- 
ulation (SSP) starting from PEGASE.2 (Fioc & Rocca- 
Volmerange (1997)) for different ages (1 Myr-20 Gyr) and 
different metallicities {Z/Zq= 0.005, 0.02, 0.2, 0.4, 1., 2.5, 
5). This library contains the sole stellar component. We 
then add the extinction based on local star-forming galax- 
ies computed by Calzetti et al. (2000) depending on Nh 
and the gas metallicity. The energetic balance allows us to 
compute the luminosity Lm, re-emitted in the infrared by 
the dust. In a similar fashion as Guidcrdoni et al. (1998) 
and Devriendt et al. (1999) , we re-distribute this en- 
ergy between 3 components (namely big grains, very small 
grains and PAH), following Desert et al. (1990). Last, 
the synchrotron emission is also included with a ly^ law 
{P = -0.7) following Condon et al. (1991). 



2.2. Free parameters 

This minimal modelling already has 6 free parameters: 

— Mssp '■ the mass of stars 

— Zssp ■ the metallicity of the stars 

— Mgas ■ the mass of gas 



- t 



the metallicity of the gas 



ssp ■ the age of the starburst 
— dsize ■ the size of the galaxy or starburst encompassing 
the previous gas and stars 
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Figure 3. Same as Fig. |^ but for a moderate starburst. The 
characteristics are taken from NGC 6090 (at z = 0.03J as 
published by Calzetti et al. ( 2000\ . Galaxies up to a redshift 



of 2 can then he studied. 



2.3. Modelling known galaxy spectra 

In order to test the validity of this simple model, we fit 
known spectra of starburst galaxies. Given the small num- 
ber of available measurements, we deliberately keep the 
metallicities to solar values. Firstly, we fit the optical part 
(< 2/im) in order to adjust M^sp, Mgas and tssp- Secondly, 
the infrared/millimetre/radio measurements are used to 
define d^ize- 

The first example corresponds to a strong starburst 
with SFR~370MQyr-i (see Fig. | and |) and the second 
one is a more moderate one with SFR~3OM0yr~^ (see 
Fig. 1^). Although real galaxies are probably composed of 
old stellar populations and young starbursts, this simple 
model based on starburst assumptions is already a good 
fit. 
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2.4. Application to FIRST 

We use the spectra of Fig. || and ^ to study the sensitiv- 
ity of PACS and SPIRE to high-z galaxies. These figures 
show starburst spectra at different redshifts as presented 
by Combes et al. (1999) together with the sensitivity of 
FIRST, ALMA and NGST. The wavelength coverage of 
PACS and SPIRE will be ideal to sample the maximum 
of the spectral luminosity distribution of starburst galax- 
ies. These instruments will be unique to understand the 
luminosity budget and properties of the galaxies repre- 
sentative of the peak of the Madau curve (Madau et al., 
1998), as well as the excitation conditions, which deter- 
mine the infrared/millimeter flux. Indeed, as shown in 
Figures |2| and ^, a galaxy at redshift 5, detected in the 
optical by the NGST, might escape detection with FIRST 
if the starburst is not strong enough. In parallel, ALMA 
will be able to detect starburst galaxies up to 2 = 30 if 
they exist, though it will not be able to sample the peak 
below z = 5. 

3. n-body simulations of galaxy mergers 

3.1. Principle 

We use these analytic starburst spectra to produce realis- 
tic maps of galaxy mergers, based on N-body simulations. 
We start from one N-body -I- hydrodynamics simulation 
(tree-code -|- SPH) of a galaxy encounter stored for dif- 
ferent epochs. It is composed of 8 000 dissipative parti- 
cles (hybrid: star or gas), 8 000 stellar particles and 8 000 
particles for dark matter. We account for star formation 
with a Schmidt law {pstar oc p^as) for the hybrid particles, 
which thus encompass a variable fraction of gas and stars. 
In order to account for supernovae, the instant recycling 
approximation for the gas rejection into the interstellar 
medium is used with standard values for the yield. 

Figure ^ displays preliminary simulations without the 
star formation and metallicity enrichment prescription for 
different epochs as well as different orientations. Although 
a single simulation can already be used to simulate a large 
variety of morphologies, we plan to generate several such 
simulations in order to span a reasonable range of possible 
morphologies. 

3.2. Application 

Figure |^ presents a very preliminary application of such 
a simulation for different epochs. Maps of such a merger 
at z = 1 are displayed as it could be observed by FIRST 
(PACS & SPIRE), NGST or ALMA. Whereas the NGST 
will be able to study the detailed morphology and possi- 
bly tidal arms, ALMA, after Ih of integration time, will 
only be able to detect the multiplicity of the source. The 
instruments of FIRST will sample the luminosity peak of 
the spectral energy distribution but will not be able to 
distinguish any signature of merging. 



We currently extend this kind of simulation in a sys- 
tematic way (different morphologies) and from a cosmo- 
logical perspective (within GalICS framework) in order 
investigate the current open questions. Firstly, the galax- 
ies that will be observed in the optical (e.g. with the 
NGST) will not be necessarily accessible to the IR/mm in- 
struments, especially if the star formation activity is mod- 
erate. Secondly, the difhculties encountered to find optical 
counterparts for SCUBA sources are hoped to be overcome 
with the next generation of instruments. This simulation 
work will help to understand the multi-wavelength depen- 
dence of these high-z galaxies and the process of galaxy 
formation. 




Figure 5. Preliminary projection of a merger at z = 1 
(LiR = 2 X 10^^ Lq) on a pixel grid with the characteris- 
tics of FIRST/PACSk SPIRE, NGST and ALMA. These 
images correspond to Ih integration time and account for 
the expected rms intrumental noise for each intrument. 



4. GalICS model 

The "hybrid" approach is a powerful tool to describe hier- 
archical galaxy formation. The GalICS model (for Galax- 
ies in Gosmological Simulations^ Hatton et al. 2001, De- 
vriendt et al. 2001) uses a large cosmological iV-body 
simulation (256^ particles in a 150 Mpc box within a 
ACDM cosmology) achieved on the CRAY-T3E of the 
IDRIS computation centre. Halos with more that 20 parti- 
cles are identified by the usual friend-of-friend algorithm, 
and a catalogue of haloes with a minimum halo mass of 
1.7 X IO^^Mq is built for each of the ~ 100 time outputs. 
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Figure 4- Sequence of a simulated merger, showing the evolution with a timestep of 100 Mgr. Each panel presents a pair 
of galaxies in interaction with 3 different orientations, as follows: upper right— plane (XY), upper left— plane (ZY), 
bottom right=plane (XZ). The green points correspond to the dark matter, the blue ones to the gas and the red ones 
to the stars. 



Merging history trees are built from these 100 catalogues. 
Finally, the evolution and merging history trees are built 
for the galaxies. The processes of dissipative gas cooling 
and collapse, star formation, chemical evolution, stellar 
feedback are modelled along lines which are detailed in 
the papers. The UV to mm spectral energy distributions 
are computed from the Stardust spectra (Devriendt et 
al. 1999). After halo merging, the satellite galaxies spiral 
down to the central galaxy with the dynamical friction 
time scale. Satellite-satellite merging is also taken into 
account. During the major mergers, a starburst light up 
and radiates mostly in the IR/submm range, and a bulge 
eventually forms. A new disk can form around this galaxy 
if it is at the cooling centre of its halo. The morphologies 
of the galaxies are fixed according to the i?-band bulge 
to disk ratios. The final output at z = includes about 
70000 galaxies, and the mass resolution progagates as an 
absolute magnitude resolution of Mb = —18. Fake images 
of deep fields can be achieved at any wavelength from the 
UV to the mm range from the output catalogues. 

So far, the code includes only very crude recipes on 
how galaxies merge and how merging triggers a starburst 
event. As a first step, it is planned to generate fake im- 
ages of deep galaxy fields at optical and IR/submm wave- 
lengths from the catalogue of morphologies of interacting 
and merging galaxies described in section 3, by choosing 
the closest case to the initial conditions for the encounter 
fixed by GalICS. As a second step, the tree-code -I- SPH 



hydrodynamical simulations can be directly run from the 
galaxy merging trees with the initial conditions of Gal- 
ICS. Such a program will allow us to assess the properties 
of mergers in a cosmological framework, and to study how 
this properties evolve with redshift. These fake fields will 
be used to study the observational strategy for deep galaxy 
surveys with PACS and SPIRE. 
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